New Limit for the Family-Number Non-conserving Decay ^ 



e T 7 



os 

OS 

os 



> 

CO 

o 
in 
o 

Os 
On 

">< 
<D 

43 



X 



M.L. Brooks, 6 Y.K. Chen, c M.D. Cooper, 6 P.S. Cooper, 6 M. Dzemidzic, 6 A. Empl, 6 C.A. Gagliardi,9 G.E. Hogan, e 
E.B. Hughes/ E.V. Hungerford III, 6 C.C.H. Jui/ J.E. Knotty D.D. Koetke,' 1 M.A. Kroupa, 6 K.A. Lan, 6 
R. Manweiler,' 1 B.W. Mayes H, 6 R.E. Mischke, 6 L.E. Piilonen,-' T.D.S. Stanislaus, 6 -' 1 K.M. Stantz, d 
J.J. Szymanski, d ' 6 R.E. Tribble, 9 X.L. Tu, s L.A. Van Ausdeln ,<? W.H. von Witsch, 6 S.C. Wright, a and K.O.H. Ziock* 

(MEGA Collaboration) 

a University of Chicago, Chicago, IL 60637 
b Fermi National Accelerator Laboratory, Batavia, IL 60510 
c University of Houston, Houston, TX 77204 
d Indiana University, Bloomington, IN 46383 
e Los Alamos National Laboratory, Los Alamos, NM 87545 
J Stanford University, Stanford, CA 94305 
9 Texas A & M University, College Station, TX 77843 
h Valparaiso University, Valparaiso, IN 46383 
1 University of Virginia, Charlottesville, VA 22901 
3 Virginia Polytechnic Institute and State University, Blacksburg, VA 24061 

(February 7, 2008) 

An experiment has been performed to search for the muon- and electron-number non-conserving 
decay fi + — > e + 7. The upper limit for the branching ratio is found to be T(n + — » e +r y)/ V(n + — » 
e+vv) < 1.2 x 10" 11 with 90% confidence. 
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It is generally believed that the standard model of elec- 
troweak interactions is a low-energy approximation to a 
more fundamental theory. Yet there is no clear experi- 
mental evidence either to guide its extension to additional 
physical processes or to predict the model parameters. 
One of these model assumptions is lepton family-number 
conservation, which has been empirically verified to high 
precision but is not a consequence of a known gauge 
theory. Indeed many theoretical extensions to the stan- 
dard model allow lepton-family-number violation within 
a range that can be tested by experiment |Q . 

The predictions of the rate for a given family-number 
non-conserving process vary among these extensions, and 
the most sensitive process depends on the model. Many 
possibilities have been explored, and the present exper- 
imental limits for a wide variety of processes have been 
tabulated in Ref. Of these, the rare muon decays 
have some of the lowest branching-ratio (BR) limits be- 
cause muons can be copiously produced and have rela- 
tively long lifetimes. The rare process, fi + — > e + 7, is 
the classic example of a reaction that would be allowed 
except for muon and electron number conservation; the 
previous limit on the branching ratio is BR(/i + — > e + j) 
< 4.9 x 10~ n j|. This decay is particularly sensitive 
to the standard model extension that involves supersym- 
metric particles [0. 

We report here a new limit for the BR of the decay 
/i + — ► e + 7 from the analysis of data taken by the MEGA 
experiment at the Los Alamos Meson Physics Facility, 
LAMPF. The dominant source of background in high- 



rate n + — > e + 7 experiments is random coincidences be- 
tween high-energy positrons from the primary decay pro- 
cess, [i + — > e + isV, and high-energy photons from internal 
bremsstrahlung (IB), fi + — » e +r yi/i7. MEGA isolates the 
fi + — » e + 7 process from the background by identifying 
the signature of the process: a 52.8-MeV photon and a 
52.8-MeV positron that are aligned back to back, in time 
coincidence, and arise from a common origin. Therefore, 
quality position, timing, and energy information are cru- 
cial. In comparison to the detector used to set the pre- 
vious limit ||, the MEGA detector sacrifices larger ac- 
ceptance and efficiency for better resolution, background 
rejection, and rate capability. It has been described in 
several papers (4j-[|] and will be discussed only briefly 
below. 

Muons for the experiment are provided by a surface 
muon beam at the stopped muon channel at LAMPF. 
The muons, which are nearly 100% polarized, are brought 
to rest in a 76 /Ltm Mylar foil, centered in the 1.5-T mag- 
netic field of a superconducting solenoid. The angle be- 
tween the muon beam and the normal to the target plane 
is 82.8° so that the stopping power in the beam direction 
is increased, while the thickness of material presented to 
the decay positrons is minimized. A sloped target plane 
also extends the stopping distribution along the beam, 
enhancing the sensitivity of the apparatus to the mea- 
surement of the decay position, which is the intersection 
of the outgoing photon and positron trajectories with the 
target foil. 

The positron and photon detectors are placed in the 
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1.8-m diameter and 2-m axial length bore of the solenoid. 
Decay positrons from stopped muons are analyzed by a 
set of high-rate, cylindrical multiwire-proportional cham- 
bers (MWPC) surrounding the target. They consist 
of seven MWPCs arranged symmetrically outside of a 
larger MWPC, coaxial with the central axis of the beam. 
These MWPCs have a thickness of 3 x 10~ 4 radia- 
tion lengths, minimizing energy loss while maintaining 
high acceptance and efficiency under the stopping rates 
of the experiment . The azimuthal location of a pass- 
ing charged particle is determined by anode wire read- 
out. The position of an event in the axial direction is 
obtained from the signal induced on stereo strips scribed 
on the inner and outer cathode foils of the MWPCs. The 
positrons come to rest at either end of the spectrometer 
in thick, high-Z material after passing through a barrel of 
87 scintillators used for timing. Outside these MWPCs, 
photons are detected in one of three coaxial, cylindrical 
pair spectrometers Each pair spectrometer consists 
of a scintillation barrel, two 250-/im Pb conversion foils 
sandwiching an MWPC, and three layers of drift cham- 
bers, with the innermost having a delay-line readout to 
determine the axial position of a hit. 

The hardware trigger, consisting of two stages of 
specially-constructed, high-speed logic circuits, is fed sig- 
nals from each of the three photon spectrometers [^J. 
Using pattern recognition programmed on the basis of 
Monte Carlo (MC) simulations, the trigger requires an 
electron-positron pair that can be potentially recon- 
structed as arising from a photon of at least 37 MeV. 
Since the instantaneous muon stopping rate in this ex- 
periment is 250 MHz, with a macroscopic duty cycle of 
6-7%, the positron chambers and scintillators have too 
many hits at any given time to be part of the trigger. 
Signals are digitized in FASTBUS with 6% dead time 
at the instantaneous trigger rate of 18 kHz. Between 
each macropulse (120 Hz) of the accelerator, the data 
are read into one of eight networked workstations, where 
an on-line algorithm reduces the data rate for storage on 
magnetic tape to roughly 60 Hz. 

Each event is characterized by 5 kinematic parame- 
ters: photon energy (i? 7 ), positron energy (E e ), relative 
time between the positron and photon (t e7 ) at the muon 
decay point, opening angle (6 e ~), and photon traceback 
angle (A9 Z ). These properties, in conjunction with the 
detector response, determine the likelihood that a signal 
is detected. The determination of the detector accep- 
tance and response functions relies on a MC simulation 
to extrapolate from experimental input to the kinematic 
region of the fi + — > e + 7 signal. To verify the MC calcula- 
tion, a number of auxiliary measurements are performed. 
The two most important are the ^ toppe ^P —* K°n — > 7771 
process and the prompt e — 7 coincidence signal from the 
IB decay. 

Pion capture at rest on hydrogen produces photons 
with energies between 54.9 and 83.0 MeV. Under the 



60 
50 
40 
1 30 
20 
10 




1 




1I1 »J*Vi*tnii 1 1 


I 

1 I7i\ U 1 1 1 1 1* 1 tJ 1 



35 40 45 50 55 60 65 70 
Photon Energy (MeV) 



FIG. 1. The E-y spectrum from photons converting in the 
outer layer of lead. The data points are produced by stopping 
pions in CH2 via the reaction iv~p — > Tv°n — > ■yjn, scaled from 
54.9 to 52.8 MeV. The curve is the response function gener- 
ated from the MC and used in the analysis of the fi + — » e +r y 
data. 



condition that the two photons have a minimum open- 
ing angle of 173.5°, these photons are restricted to have 
energies close to 54.9 and 83.0 MeV and a spread much 
smaller than the detector response. Figure [l] shows the 
experimental line shape for the 54.9 MeV photon for con- 
versions in the outer Pb foils of the three pair spectrome- 
ters, scaled to 52.8 MeV. The curve is the response func- 
tion generated from the MC that is used in the analy- 
sis of the /i + — > e + 7 data. We attribute differences in 
the low-energy tail to charge exchange of in-flight pions 
from carbon in the CH2 target and discrepancies in the 
high-energy tail to contributions from other opening an- 
gles due to reconstruction problems for the high-energy 
photon. The measured and simulated line shapes agree 
better for conversions in inner Pb foils, which have worse 
resolution. The energy resolutions are 3.3% and 5.7% 
(FWHM) at 52.8 MeV for conversions in the outer and 
the inner Pb layers, respectively. The tt° decays also pro- 
vide the time response between the two photons, which is 
reasonably characterized by a Gaussian with a a = 0.57 
ns for each photon. 

Observation of the IB process demonstrates that the 
apparatus can detect coincident e — 7 events. At nomi- 
nal beam intensity, this process is completely engulfed by 
random coincidences. Figure shows the spectrum for 
i e7 , with the beam intensity reduced by a factor of 60, the 
magnetic field lowered by 25%, and the fi + — ► e+7 on- 
line filter suppressed. The peak shown is for all energies 
of the detected decay products. The area of the peak is 
very sensitive to the exact acceptances of the detector at 
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FIG. 2. Values for t el from the process — » e + ^vV under 
the conditions of reduced rate and magnetic field. 

its thresholds and can be calculated by MC simulation 
to within a factor of two. If the data and the simula- 
tion are restricted to E 1 > 46 MeV, E e >40 MeV, and 
6 ei > 120°, the branching ratio is reproduced within un- 
certainties. The shape of the peak can be characterized 
by a Gaussian with a a = 0.77 ns. The dominant contrib- 
utor is the photon timing, as measured in the stopping- 
pion experiment, which must be scaled down from about 
70 to 40 MeV for the comparison. At 52.8 MeV, the MC 
simulation indicates the photon-positron resolution is a 
= 0.68 ns. 

In the IB and fi + —> e + 7 processes, the origin of the 
photon is defined to be the intersection of the positron 
with the target. The photon traceback angle, A9 Z , speci- 
fies the difference between the polar angles of the photon 
as determined from the line connecting the decay point to 
the photon conversion point and from the reconstructed 
e + — e~ pair. The resolution of A6 Z is dominated by 
multiple scattering of the pair in the Pb converters. The 
observed response for inner and outer conversion layers 
for the IB process is in excellent agreement with the MC 
simulation. The traceback resolutions appropriate for the 
/i + — > e + 7 analysis are a — 0.067 and 0.116 rad for con- 
versions in the outer and the inner Pb layers, respectively. 

The resolution of E e is determined by the slope of the 
high-energy cut-off edge in the spectrum of the decay, 
[i + — > e + uV. It depends on the "topology" of the track, 
which is determined by the number of loops these parti- 
cles make in the magnetic field between the target and 
scintillator and the number of chambers they traverse. 
The E e spectrum is shown in Fig. || for one of three topol- 
ogy groups. The MC line shape is characterized near the 
centroid by a Gaussian and in the tails by different pow- 
ers of the deviation from the central energy. To extract 



the response function from the data, this line shape is 
convoluted with the spectrum from normal muon decay, 
modified by detector acceptance and unphysical "ghost" 
tracks. Ghost tracks are a high-rate phenomenon and are 
reconstructions made from the fragments of several phys- 
ical tracks. They are the source of events well above the 
kinematic limit for the positron energy. The solid curve 
in Fig. H is the fit, and the dashed curve is the corre- 
sponding line shape. The central Gaussians of the three 
topology groups have a = 0.21, 0.23, and 0.36 MeV. 

There is no way to measure the response function for 
6l e7 . The MC simulation is relied upon to produce this 
distribution and gives the FWHM for cos(0 e7 ) as 1.21 x 
10~ 4 at 180°. Given helical tracks, knowing the location 
of the target is critical to obtaining the correct absolute 
value of 9 ej , and the mechanical survey provides the most 
accurate measurement for the analysis. 

The data for this experiment have been taken in three 
calendar years, 1993-95. The full data set is based on 
1.2 x 10 14 muon stops collected over 8 x 10 6 s of live time 
and results in 4.5 x 10 8 events on magnetic tape. These 
events are passed through a set of computer programs 
that reconstruct as many as the pattern recognition al- 
gorithms can interpret. The programs include physical 
effects such as mean energy loss in matter and non- 
uniformities in the magnetic field. Events are required 
to satisfy separate xt cu ^s on the positron and photon 
fits and loose cuts on the signal kinematics (E e > 50 
MeV, E 1 >46 MeV, \t n \ <4 ns, cos((9 e7 ) < -0.9962, and 
|A0 Z | < 0.5 rad). Events in which the positron momen- 
tum vector at the decay point appears to lie within 5° of 
the plane of the target are discarded. After roughly one 
year of computing on a farm of UNIX workstations, the 
data set has been reduced to 3971 events that are fully 
reconstructed and of continuing interest. This sample is 
large enough to allow a study of the background. To re- 
move incorrectly reconstructed events, the images of the 
photon showers in the pair spectrometers are manually 
scanned. The efficiency for keeping real photons is moni- 
tored by mixing about 500 52.8-MeV MC events into the 
sample in a non-identifiable way and finding that 91% of 
the MC events pass, whereas only 73% of the data events 
are selected. 

The acceptance of the apparatus is obtained by simu- 
lating 1.2 x 10 7 unpolarized /i + — > e + 7 decays and find- 
ing that 5.2 x 10 events survive processing by the same 
codes used for the data analysis. Thus the probability 
that a /i + — > e + 7 decay would be detected is 4.3 x 10~ 3 . 
This value is reduced by 20% to account for inadequacies 
in the MC simulation that over estimate the acceptance. 
The shortcomings primarily involve inter-channel cross 
talk and are estimated by comparing the images of many 
data and MC events. The acceptance is further reduced 
by 9% for the inefficiency of manual scanning. The total 
number of muon stops is determined by calibrating the 
rates in the positron scintillators to a known muon flux. 
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FIG. 3. The E e spectrum from ji + — » e + vV extracted from 
full rate data for the middle topology group. The solid curve 
is the fit used to extract the line shape (dashed curve). 



After correcting for dead time, the single event sensitiv- 
ity for the experiment is 2.3 x 1CP 12 = l/A^, where JV M 
is the number of useful stopped muons. 

The determination of the number of fi + — ► e + 7 events 
in the sample is evaluated using the likelihood method 
described in the analysis of previous experiments . The 
formula for the normalized likelihood is 
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where A — 3971, A e7 is the number of signal events, Nib 
is the number of IB events, and P, Q, and R are the prob- 
ability density functions (PDF) for signal, IB, and ran- 
doms of each of the five parameters describing the event. 
The PDFs P and R are the products of statistically inde- 
pendent PDFs for the five parameters, each normalized 
to unit probability over the full range of the variable for 
the sample. The signal distributions are taken from MC 
distributions as described. The background PDFs are ex- 
tracted from the spectral shapes of a much larger sample 
of events, where the constraints on the other statistically 
independent parameters remain very loose. Here Q is 
taken from MC simulation of the IB and has correlations 
amongst the variables. The events fall into the following 
categories: positron topology, photon conversion plane, 
target intersection angle, and data taking period. As a 
result, PDFs are extracted for each class of events and ap- 
plied according to the classification of individual events. 

The likelihood function evaluates the statistical separa- 
tion between signal, IB, and background. To observe the 
impact of quality constraints in the pattern recognition, 
they have been relaxed to produce a sample three-times 



larger. One event emerges with a large value of Pj R that 
is significantly separated from the distribution. However, 
this event has a large positron xti indicative of a ghost 
track. The adopted constraints produce a sample with 
considerably less background. The peak of the likelihood 
function is at A e7 =0 and Ajb=30±8±15. The system- 
atic error assigned to Nj b is due to the uncertainty in the 
shape of the background time spectrum when the events 
arc filtered by the online program. The expected number 
of IB events is 36±3±10, where the systematic error is 
due to finite resolution effects across the cut boundaries. 
The 90% confidence limit is the value for A e7 where 90% 
of the area of the likelihood curve lies below A e7 and 



Nib is maximal. This value is A e7 
limit on the branching ratio is 



< 5.1. Therefore, the 
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In comparison to the previous experimental limit || , this 
result represents a factor of 4.1 improvement. The pre- 
vious experiment would have had 58 background events 
at the same sensitivity instead of the 2 found here. This 
improvement further constrains attempts to build exten- 
sions to the standard model 0. 
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